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We study the effect of a single excluded site on the diffusion of a particle undergoing random 
walk in a d-dimensional lattice. The determination of the characteristic function allows to find 
explicitly the asymptotical behaviour of physical quantities such as the particle average position 
(drift) {x){t) and the mean square deviation (x^)(f) — (i)^(f). Contrarily to the one-dimensional 
case, where {x){t) diverges at infinite times {{x){t) ~ f^") and where the diffusion constant D is 
changed due to the impurity, the effects of the latter are shown to be much less important in higher 
^r^ ' dimensions: for d > 2, {x){t) is simply shifted by a constant and the diffusion constant remains 

^^ I unaltered although dynamical corrections (logarithmic for d = 2) still occur. Finally, the continuum 

space version of the model is analyzed; it is shown that d = f , is the lower dimensionality above 
which all the effects of the forbidden site are irrelevant. 

Oh! 
(D . 

C/^ ' I. INTRODUCTION 

^~^ ' InteractioHS between diffusing particles are expected to play an important-, role in many systems of physical 

, ■ interest (zeoliteail, biological membrancsofl, one-dimensional hopping conductivityo) , especially in low dimension and 

►^ [ even with short-range interactions. For d = l|-|the case of a contact interaction has been extensively studied in the 

ly-v i past, following the pioneering paper by HarrisO on the so-called tracer problem, first showing that the mean-square 

f — displacement (a;^) increaseSras t^'^ at large times. This result was also obtained by van Beijeren et alB and, more 

O^ . recently, by Rodenbeckei aln; these authors indeed solve the full A^-particle problem for an arbitrary initial condition 

0^ ■ and recover Harris' result by going to the TV —> -f cx) limit. Still in the case of a zero-range repulsive interaction, 

^^ , asymptotic results (transport coefhcients and distribution laws) in the extreme case of a finite- A^ compact initial 

^^ '■ cluster have been obtained in refcl, where the relation with the theory of extreme events was also discussedB. 

— .^ ] Obviously, the effects of short-range interactions are expected to become less and less important in higher 

jrt ' dimensions (on the other hand, the interplay between interaction range and space dimensionality remains, as far as 

we know, an open question). In the present paper, we address the simplest problem, namely that of two particles 

with a contact interaction undergoing d-dimensional lattice random walk. After transforming to the center-of-mass 

'^ ' frame, this essentially reduces to the Brownian motion of a single particle subjected to an isolated reflecting barrier. 

^ The equations for the general d-dimensional case are given below and, among other things, allow to show that d = 1 

O ' is indeed the marginal dimensionality above which the random walk is unaffected by the presence of the localized 

impurity in the continuous version of the model; by this, it is meant that transport coefficients, as deduced from the 

dominant term in the asymptotic regime of physical quantities, are the same with and without impurity for d > 1. 

Nevertheless, corrections still occur in the lattice model; that these corrections, less and less relevant as d increases, 

basically display a a'^ behaviour for d > 1, where a denotes the lattice spacing. 

C^ ' This paper is organized as follows. We first write down the basic equations for any dimension d and give 

the expression of the Laplace transform of the generating function giving by derivations all the moments of the 

particle coordinate in the lattice. We then analyze this general result, retrieving well-known results for the d = 1 

case. Subsequently, we examine in details the two-dimensional case, which exhibits lattice-dependent logarithmic 

corrections, and give the main results for the d > 2 case. Even in such high dimensions, the presence of the forbidden 

site has consequences on the asymptotical dynamics: the final average particle position displays a shift Sx{ which 

scales as a'^x^'^'^^ where xq is the distance between the starting point of the particle and the exclusion site. Finally, 

the space-continuous version solution is written in full, showing that the moment-generating function is, in this limit, 

insensitive to the presence of the excluded site as far as the dimensionality d is strictly greater that one. 
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II. BASIC EQUATIONS FOR THE D-DIMENSIONAL LATTICE WALK 

In the following, we consider a particle undergoing random walk in the viscous limit on a d-dimensional hyper- 
cubic lattice having a forbidden site. Denoting Cp the unit-vector in the p-th direction of the lattice, an arbitrary site 
of the lattice is fully specified by its vector n = X)p=i ^p^p! the excluded site is located at the origin. Let P{n, t) be 
the probability to find the particule at site n at time t. In the continuous time description, it is readily seen that the 
time evolution of P{n, t) is governed by the following master equation: 



— P(n, t) = W[l - d{n, 0)] J2 "^ ^^ ~ ^^"' + ^^P' ^)] [P(n + eep, t) - P{n, t)] 



di 



(1) 



£=±lp=l 



where W denotes the probability per unit time to jump from one site to any of its nearest-neighbours; in (nl), 5{n, m) 
equals 1 if n = m and vanishes otherwise; obviously, P{n = 0,t) = at all times. The easiest way to solve this 
problem is to deal with the characteristic function of the probability distribution: 



i:{$,t) =^e'*~-"P(n,t) , 



(2) 



which allows to find the probability function P by the inverse formula : 



(3) 



where (jjp = 4>-(ip and where / A'^'cj) means {2tt) '^ J_^ d(l)i J_^ d(j)2 ■ ■ ■ /_ J d^^. Multiplying both sides of (|]) by e"^-" 
and summing over n, one readily obtains: 



- V(0,t) = 2M^^(cos0p - l)V(<?,t) + W^ E E (e"*" - 1) / 'i''^' (^""^'' - l) ^('?''*) 

"^ p = l £ = ±1 p=l J ^ ^ 



(4) 



For sake of simplicity, the initial position of the walker is chosen to be xq = riQaei. With this condition, and since 
there is no external bias, the probability distribution is invariant under a mirror transformation through any plane 
perpendicular to Cp with p ^ 1. In addition, the quantities / d'^(^' cos 4>pip {4> ' ,t) , p ^ 1, are all equal. Introducing 
the Laplace transform of V'(0j t) ■ 



4'i(f>,z) 



die-^V(0,i) 



(5) 



one readily obtains: 



(d \ d 

z + 2VFE(1-cos0p)Ua(0,z) = Vo + 2VFE(cos0p-1) / d''^' (cos 0^ - 1)^^(0', z) 
p=l / p=2 •' 

+ 2W d'^cj) ' [(cos 01-1) cos 01 + sin 0i sin 0'^ -0(0 ', z) , 



(6) 

where -00 = V'(0! 0) = e'"""^^. Equation (0) is a homogeneous separable Fredholm system, and can be solved by 
quadratures (see Appendix H), yielding the expression of V'('/'; z) (see eq. ( A8). For further reference, we explicitly 
write down the Laplace transforms of the average coordinate and mean-square displacement: 



{xp){z) = -ia 



9-0(0, z) 



{xl){z) ^ -a 



2 520;(0, z) 



A=0 



(7) 



A=0 



With the solution given by (A8), one precisely has 



{xp){z)^5pi- no-\Si{Z) 



(8) 



{x^)iz) 



2 r 



+ nl + CiiZ) + id-l)C{Z) 



(9) 



where Z = z/{2W) and where the functions 5*1, Ci and C are defined in Appendix ^. In the next sections, we 
shall find the large-time behaviour of these quantities according to the dimensionality d, by analyzing their small Z 
expansions (Z < 1 ^=^ t > W~^). 



III. THE ONE-DIMENSIONAL CASE 



The reduced problem is here simply the 1 d diffusion of a particle submitted to a refiecting barrier. Although the 
result is well known, we shall discus s it t o stress the main difference between this case and that of larger dimensionality. 
For d = 1 the relev ant Ci function ( A4 ) greatly simplifies (Ci = /3„o /(l — ai — /3i ) iy,whereas Si is formally unchanged 
and still given by (A4). Let first us discuss the drift term, given by dq). By usin^: 



dx e-°"'ln{x) = 



(a + Va^ ~ l)"\/c 



axx-^e-""" 



In (a;) 



(a — y/a^ 



1 



(10) 



we find that 7„ (introduced in (A7)) is well-defined for all z, '3iz > and indeed has a finite limiting value, namely 
jn{z = 0) = 1; this is one peculiarity of the one-dimensional case, for which jn — 1 ^ V2Z at small Z, so that Si 
diverges in the limit Z — > 0. In higher dimensions, the modulus of 7„ is always stricty smaller than 1 for all Z , so 
that 5*1 no more diverges in the limit Z —^ 0. This property is crucial since it drastically changes the behaviour of 
the drift term according to the dimensionality, setting a marked difference between the cases d = 1 and d > 1 (see 
also section M, where the continuum limit is analyzed). This fact points out d = 1 as being the marginal dimension 
for the drift term. A detailed calculation yields: 



Siiz) 



Similarly: 



Ci{z) 



1 1 


-"0 + 


'no, 
y(no- 


-1)- 


1" 

" 8_ 


V2Z 2 


1 1 


-no + 


no, 
y("o- 


'!)- 


3' 

' 8_ 


V2Z ' 2 



^2Z 



'2Z 



(11) 



(12) 



Using (O) and (E4), the small-z expansions of (xi){z) and {x\){z) can now be written out. Performing then the 
Laplace inversion yields the following asymptotic expansions: 



(xi)(i) - 2a 



a {2^o_^^ I 

2 %a^nmi 



(13) 



(a;?)(i) - 2c?Wt -\- 2a^\\ V {x\ - axo + ^) + ^(4^2 - 4axo - ic?) 



2 1 

8* 



1 



yfiml 



+ o(t-i) 



(14) 



Note that the initial condition xq does not appear in the two first-dominant terms. The mean square deviation for 
the coordinate xi results: 



^x\(i)^ {x\)(i)-{{xi){i)f 



2 1 



a^Wt 



K'-'-i)' 



O(t-i) 



(15) 



The continuum limit is obtained from the lattice model by taking the limit T4^^cxD,a^0, tiq— > +00 with a}W = D 
and noa — xq finite; performing this, one finds: 



iximr.2^+^^^+0{t-') , Axlit) -2(^1-^) (^Dt+f^+Oit~') 



(16) 



Note that the continuous limit strongly modifies the asymptotic expansions: in this limit, the first sub-dominant term 
of the lattice model drops out and the first relative corrections to the asymptotic leading term are 0{t~^) instead of 



0{t~^^^) in the lattice version. On another point of view, comparison with (|lj) and (13) shows that the first lattice 
corrections vanish as the first power of the lattice constant a. 

The above expressions, obtained as the continuous limits of the lattice model, can be compared to the well- 
known results directly obtained in the continuous framework. For a particle subjected to a reflecting barrier and 
starting a distance 2:0 > apart, an elementary calculation yields the moment-generating function in the continuous 
version, which reads: 



ij{k,t) = i[V>+(fc,t)-hV-(fc,t)] i^±{k,t) 



1 + $ ±- 



^— Dk t ± ifcxo 



where $ is the Error functiont3. By successive derivations at /c = 0, the two first moments result: 



{xi){t) = Xo$ 



ADt 



-xi/(iDt) 



{xi){t) = 2Dt + x^o 



'ADt/ V TT 
These expressions are exact for all i; by an expansion at large times, they precisely reproduce (|16|), as it must. 



(17) 



(18) 



IV. THE CASE D > 2 



The expected irrelevance of the forbidden site as dimension is increased will now be explicitly displayed. Among 



other things, it will be seen that, for d — 2, 



and Ax contain logarithmic corrections to the "bare" diffusive 



regime, which become constant in time for d > 2. These results allow to state that d = 2 is indeed the marginal 
dimension for the diffusion constant characterizing the dominant term of the asymptotic behaviour. On the other 
hand, a permanent shift of the coordinate persists for any d, even at very large times. 

Let us begin with the average coordinate (xi). Contrarily to the case d = 1, the integral 71 is always strictly 
smaller than one, since e~^'^+^-'^/„(a;) < 1 for all Z and x real. Using results given in the Appendix B, one eventually 
obtains the following expression for the Laplace transform of the average coordinate valid for \z\ ^ W: 



which gives: 



{xi){i 



xq A+{2/nl)-{no-l)ZlnZ 
z 4 — no Z\nZ 



i^iKt) 



Xo 



a 

2x0 



1 - 



1 



16Wt 



xq + 6x{t) 



(19) 



(20) 




FIG. 1. Comparison for the final shift coordinate 5xf between analytical (solid curve) and numerical results (dashed line), 
for 710 = 1 and d — 2. The unit of time is W~^. 



The final value of the average coordinate is thus simply shifted by Sxf = a? /{2xq), all the more since the latter starts 
close to the forbidden site. As an illustration we show on figure 1 the analytical (solid curve) result for the final shift, 
comp ared w ith the numerical calculation (dashed curve), in the case rig — 1. The second moment can be obtained 
from ( A4), (A5) and (^, Using the results of Appendix B: 



{x^){z) 



Aa^W 



(Z + 2)a„„(Z, 2) 



l-Zai(Z, 2) 



(21) 



{x^){z) ^ - 



Aa^W 



In- 



2W 



Laplace inversion gives the dominant behaviour at large times: 

{x'^){t) r^Aa^Wt + xl + —\n2Wt 



From ( poD and (23), one finds the mean square deviation: 

Ax'^it) ^ Aa^Wt 



— In 2Wt 

n 



At this point, we see here that the diffusion constant, defined as usually: 



D 



lim 



i^^^^W ' 



(22) 



(23) 



(24) 



(25) 



is now unchanged by the impurity but a logarithmic correction occurs in the sub-dominant term, which can be viewed 
as the reminder of the change of D when d = 1: the diffusion constant is unchanged but the first correction diverges 
infinitely slowly at infinite times. This is clearly a marginal correction: at times large enough, the two-dimensional 
diffusive regime is the same with or without impurity. 

In the continuous limit, the expressions (pH) and (p4) give: 



(Xl)(i) r^ XO+O 



Ax^{t) - ADt + o(a^ In 



2a2 



(26) 



Let us now turn to the case d > 2, for which all the integrals occurring in the subdominant terms are strictly 
convergent for all Z, 5ftZ > 0, whereas the dominant terms are exactly the same as in the absence of a localized 
impurity (indeed, all the integrands behave at most as x~'^''^ at infinity). More precisely, one finds: 



and: 



{xi){z) 



{x^){z) 



xo_ 

z 



1 



C("o, d) 



l-e(l, d) 



1 l' 2da'^W 

z \ z 



A 



where A is a constant and where the function ^{n, d) is defined in (B13). It results: 



(xi)(i) 



Xq 



1 



C("-o, d) 

i-e(i,d) 



Af2(i) ^ 2da^Wt + A' 



(27) 



(28) 



(29) 



where A' is another constant. The correction to unity in the brackets clearly represents the microscopic exclusion 
phenomemon, which agai n gives a final shift of the particle position, due to the presence of the impurity. As shown 
in the Appendix B, (B17), the function ^{uq, d) is approximately given by: 



iin,d)^i^n)-^r{d/2) 



(n > 1) , (30) 

where F denotes the Euler Gamma function. Equations (|2j) and (pfl) show that the final coordinate shift has an 

algebraic behaviour (Sx ^ a** xo ) as a function of the distance between the starting point and the forbidden 

site. On the other hand, the subdominant term in Ax'^ turns out to be now a constant in time, instead of diverging 
logarithmically as is the case for d = 2. 



V. THE SPACE CONTINUOUS LIMIT 

Here and there in the above, we gave the continuum hniits of the first moments, in order to emphasize the 
corrections in the asymptotic regimes which are due to the lattice. In this section, we exphcitly show that the d > 1- 
random walk in continuum space is, as a whole, unaffected at all times by the local impurity. The central quantities 
Si and Cp only depend on the reduced variable Z = za^ /{2D), pointing out that the continuum limit embodies the 
limit Z ^r Q. Let us first consider the behaviour of Si, using the results of the Appendix B, we find at the lowest 
order in Z : 



l-7r-2r(d/2) -Fd{Z) 



where: 



Fi(Z) ~ -V2Z , F2{Z) ~ ^Z\ogZ , Fd>2{Z) ~ Cte x Z . (32) 

Here we clearly see that for the drift term, the marginal dimension is d = 1: the term Si is divergent at low Z only 
for d = 1 and remains fi nite in this limit for d > 1. Since the drift only depends on Si, we can conclude that, for 
d = 1 , {xi){t) — xq = 2^J Dt/iT, whereas for higher dimension, {xi) — xq vanish as a in the cont inuu m limit. 



Generally speaking, it is easy to obtain the continuous limit of the characteristic function (A8). A few algebra 



yields 

^{k, z) = ^ ^^^^ (V^o + ifc\/ 7 + g(«)) [d^l), (33) 

^k, z) = JT^CV-o + 0(a'^)) (d > 1) . (34) 

Thus, when the continuum limit is performed, all the additional terms arising from the excluded site strictly vanish 
except for d = 1. As a whole, it turns out that d = 1 is the marginal dimension of the problem in its space continuous 
version. 

As a final remark, let us note that the present work indeed solves the problem of two particles with a hard-core 
mutual repulsion (exclusion process), each of them having the same diffusion constant D. By separating the free 
diffusion of the center-of-mass - which undergoes a normal diffusion with a constant D jl -, one is left with a reduced 
particle with a constant 2D and subjected to a static excluded point, which is the problem fully solved above. 

VI. CONCLUSIONS 

We have studied the effect of a forbidden site on the random walk on a hypcrcubic lattice according to the 
dimensionality d of the latter. For the mean square dispersion of the coordinate, the diffusion constant D defining 
the asymptotical dynamics is altered only when d = 1; for d = 2, logarithmic corrections ~ Ini occur, which become 
constant in time when d > 2, whereas D is unchanged. On the contrary, the average coordinate is sensitive to the 
impurity in all cases. For d = 1, it goes to infinity at large times {{x){t) ^ t^^^)', for d > 2, its limiting value is 
always finite and displays a shift 5x{ ex u'^Xq '^^ as compared to its initial value xq. 

As expected, these effects are found to be less important in the continuous- space version of the problem. By 
considering the generating function of the moments, it was shown that the impurity is totally irrelevant except for 
d = 1: as far as d > 1, the continuous-space generating function is exactly the same as for the ordinary random walk 
(the above- mentionned corrections for d > 2 are thus lattice effects). This allows to state that d = 1 is the marginal 
dimension for the continuous-space problem. 

All these results were obtained by assuming a contact interaction between the walker and the impurity. Obvi- 
ously enough, it can be anticipated that a long-range interaction strongly alters the present conclusions. Up to our 
knowledge, the interplay between the dimensionality and the interaction range is, presently, an open question that 
clearly deserves further study. 

We are indebted to Julien Vidal for his most valuable remarks and comments on the manuscript. 



APPENDIX A: SOLUTION OF THE FREDHOLM EQUATION 

Equation (^) is a homogeneous separable Fredholm system and can be solved by quadratures. We set: 

z + 2W}2p=ii^- COS (l)p) 
and introduce the auxiliary dimensionless quantities: 

Cp{z) = 2W I d'^ (t) cos (j)p'4!{$,z) , Si{z) = 2W [ d'^cj) sm(t)iiP{(f,z) . (A2) 



Due to the chosen initial condition, the (rf — 1) functions Cp, p = 2, 3, . . . , d, are all equal and arc simply denoted 
by C in the following. Using the standard procedure for solving such a separable Fredholm equation, Cp and Si are 
seen to be given by an inhomogeneous system, which can be readily written and solved. Then, using the identity 
u~^ = L dxe""^, 3fJw > 0, it turns out to be that all these quantities can be expressed in terms of integrals 
involving the Bessel functions of the second kind /„. More precisely: 

an{Z,d)= dxe-(^+'^)^I„(x)lg-l(x) bnn'{Z,d)= dx e-^^+''^^ In{x) In' {x) I^^'^x) , (A3) 

JO Jo 

A somewhat lengthy but straightforward calculation gives: 

(7i - ^ [1 + ai - /3i - (d - 2)(an - ai)]/3„„ + {d - l)(an - ai)a„„i] , Si = -^^^ , (A4) 

D 1 - 7i 

Cp^i = C = —[{1 + ai- /3i)a„oi + (an - ai)/3„J . (A5) 

The various quantities appearing in these expressions are defined as follows: 

D = [l + ai-(3i-{d- 2){aii - ai)](l + ai - Pi) - {d - l)(aii - aif] , (A6) 

Pn = ^[""-1+^n+l] ' 7n = l^ [a-n-l - ttn+l] , «„„' = 6„„'(Z, d-1) . (A7) 

With these definitions, the solution of the central equation (^) writes: 

d 



-ipo + C{z) ^(cos(/ip - 1) +Ci{z) (cos0i - 1) + 5*1(2) sin 01 

p=2 



i;i^,z) = G(0,z) 
and is the starting point of the analysis undertaken in the present work. 



(A8) 



APPENDIX B: ASYMPTOTICAL BEHAVIOUR OF MANY-BESSEL FUNCTIONS INTEGRALS 

In this paper we repeatedly have to estimate integrals of the kind: 

/•OO 

I{Z,d,n,m)= dxc-^''+^>InAx)...InJx) x^ , (Bl) 



in the limit Z ^ 1, where In{x) denotes a modified Bessel function and n — {ni, . . . , Ud}- When Z = 0, no 
singularity arises from the lower bound provided that ni + . . . + n^ + to > —1 ( this is always the case in this paper). 
On the other side {x -^ 00), the integral is finite for m < {d/2) — 1. Our goal is to obtain the asymptotical behaviour 
of /(Z, d, n, m) for large rii (but numerical checks show that some of the following analytical results still hold true 
even for n? — Sti^ = !)• 

In particular, when I(Z,d,n,m) is divergent, the knowledge of the nature of the divergence for Z ~ is 
required. Let us first introduce the generalized integrals : 



Ig{Z, d, n, m,Lj)= dx c-(^+'^)"e-"/"/„,(x).../„,(x):E™ . (B2) 

Jo 

which share the same diverging properties than / when Z ^ 1 and allow to find the latter by one or other of the 
following procedures: 



I{Z,d,n, m) = lim Ig{Z, d, n, m, to) 



(B3) 



Ig{Z, d, n, m— \, uj) 



d_ 



Ig[Z, d, n, m, w) 



(B4) 



Ig{Z, d, n, m + 1, Lu) = -Tr^/g(^, d, n, to, to) . 



d_ 
dZ' 



(B5) 



Let us first examine the case d even {d = 2p). The first non- logarithmically diverging integral is Ig{Z, 2p, n^p, oj). 
The divergence near Z = arises from large values of x. In order to get for purpose a convenient asymptotical form of 
Bessel functions, we first use an integral representation of the latter, and approximate it by using a stationary-phase 
argument when both n and x are much greater than unity : 



1 r 1 r°^ 

27r J^rr 27r J_^ 



1 



/27ra; 



(x, n > 1) 



(B6) 



This approximation being done, all the integrals Ig and / onl y depend on the length p£ the vector n, simply denoted 
by n in the following. Using the approximate expression (B6), we obtain with help oiIj : 



Ig{Z, 2p, fi, p, u) ~ (27r)~ 



dx 



-Zx^ 



o + n^/2 



{2.)-^^l^±f^K,{vw^T^mz) , 



(B7) 



Ki denoting the Bessel function of the second kind. 

For d odd {d — 2p + 1) we choose Ig{Z, 2p + 1, n,p, uj) as a generator; by the same way as above we find: 



Ig{Z, 2p+l, n, p, cj)~(27r)-(P+^) 



dx e~^=" e~ 



o + n^/2 



^ = {2TT)-'-P-i\ - e-V^('^+"V2) 



(B8) 



Expressions (B7) and (B8) are correct in the limit Z <C 1 and n^ ^ 1. Using then (B4) and (B5), we are now able 
to calculate the asymptotical behaviour of the whole family of / integrals. 

As an example, let us consider I{Z, 2, {0,n}, —1), which occurs in the calculation of the average coordinate 
for d = 2. Explicitly, one has: 



/•oo 

/(Z, 2, {0,7i}, -1) = / dxx-ie-(^+2)^/o(x)/„(x) , 
Jo 

and the method sketched above yields: 

I{Z, 2, {0,n}, -1) ~ — - + —ZlnZ (Z < 1 , n > 1) 



Trn^ 2tt 



(B9) 



(BIO) 



We numerically checked the value of the Z-independent constant; it turns out that for n > 5, l/{nn'^) gives the 
correct value with a relative error smaller than ~ 10~^. On the other hand, we found no simple way to numerically 
check the validity of the first-order term in the expansion (BIC), but aae can analytically prove that it is correct for 
I{Z, 2, {0, 1}, -1). Indeed, the integral I{Z, 2, {0, 1}, 1) is computed inll3: 



/(Z, 2, {0,1}, 1) 



27r 



Z 



'2-1 ^Z'' Z' ^ Z' ' 



(Bll) 



where K and E are the complete elliptic functions of the first and the second kind respectively, and where Z' = l + Z/2. 
After a few algebra, we find : 



lim ZJ(Z, 2, {0,1}, 1) = ;^ 
z ^0 Zn 



(B12) 



Setting: 



(B12) entails: 



£,{n, d) = /(O, d, n, -1) , (B13) 



lim /(Z, 2, {0,1}, -l)=e({0,l}, 2) + -!-ZlnZ + CiZ + ... , (B14) 

z ^0 Ztt 



which for n = 1 is in agreement with (BIO). 

As an apphcation, let us find the asymptotical behaviour of the function ^{n, d) defined by ( B13| ), which is 
strongly related to the coordinate shift in any dimension. For d = 2p, we find (remember that n = \fi\J: 



e(n,2p) = (-1)"+^ 7n;TT^siZ, n, 2p, p, lu) 



dujP+ 
whereas for d — 2p + \, one has: 

C(n,2p+ 1) = (-If+i -d-i;+i^9{Z, n, 2p+l, p, u;) 



~ 7r-P(p-l)!n-2p , (B15) 

Z=ui=0 



^^_p 1.3.5... (2p-l) ^_(,p^i) ^^^^^ 

Z=uj=0 2^ 



These two definitions turn in fact into only one : 

^(n,d) ~ n-^T{d/2)n-'^ , (B17) 

r still denoting the Euler Gamma function. It is worth noting that this asymptotical behaviour, a priori only valid 
for \n\ 3> 1, is very rapidly convergent and actually holds true even for rather small values of |n |, as shown by a few 
numerical checks. 
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